Effect of a polymer additive on heat transport in turbulent Rayleigh-Benard 
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Measurements of heat transport, as expressed by the Nusselt number Nu, are reported for turbu- 
lent Rayleigh-Benard convection of water containing up to 120 ppm by weight of poly- [ethylene oxide] 
with a molecular weight of 4 x 10 6 g/mole. Over the Rayleigh number range 5 x 10 9 <, Ra <, 7 x 10 10 
Nu is smaller than it is for pure water by up to 10%. 

PACS numbers: 47.27.-i, 47.27.tc 



It has been known for about six decades that minute 
quantities of a polymer dissolved in a fluid can reduce sig- 
nificantly the drag experienced by turbulent flows near 
solid surfaces (see, for instance, [JQ). During this time 
literally thousands of papers have been written about 
this technologically important problem. However, we are 
not aware of any experimental study, and know of only 
a single very recent theoretical investigation [3[, of the 
influence of polymer additives on turbulent convection 
in a fluid confined between horizontal parallel plates and 
heated from below (known as Rayleigh-Benard convec- 
tion or RBC). In this letter we present data showing 
that the addition of about 100 parts per million (ppm) 
by weight of poly- [ethylene oxide] (PEO) can reduce the 
heat transport in turbulent RBC by 10% or more. 

The reason why the influence of dilute polymers on 
RBC has not been studied before is easy to see. In typi- 
cal flow geometries, such as pipe flow, it was found [4| 
that significant drag reduction occurs when the poly- 
mer relaxation times are comparable to or longer than 
the time scales of some of the fluctuations in the tur- 
bulent flow. Turbulent RBC (see, for instance, [1, 
generally involves only modest Reynolds numbers. Thus, 
even though there is a continuum of time scales and eddy 
sizes, the natural time scales of the most abundant ed- 
dies or fluctuations are typically of order a second or 
longer, [3, [gj, |9j and fluctuation times comparable to typi- 
cal high molecular weight polymer relaxation times (see, 
for instance, [l(|) are expected to be virtually absent. In 
addition, the viscous boundary layers (BLs) above the 
bottom and below the top plate, although fluctuating, 
are laminar because the shear Reynolds number, defined 
in terms of the BL thickness, usually is only of order 10 2 
or less. Thus there was no obvious reason to study the 
influence of minute quantities of polymers in this system. 

In view of the above discussion of the time scales of 
the problem it seems likely that the mechanism for the 
experimentally observed reduction of the heat transport 
in RBC is quite different from the usual drag-reduction 
mechanism. Since the heat transport at large Rayleigh 
numbers is dominated by the emission of plumes from the 
thermal boundary layers, the experiments suggest that 



the emission is reduced by the polymer additive. In view 
of the fact that the thermal BLs are, roughly speaking, 
marginally stable, and that plume emission may be re- 
garded as a manifestation of this near-instability, it is 
not unreasonable that this process should be exception- 
ally sensitive to such factors as polymer concentration; 
but the precise mechanism involved in this suggested 
"smoothing" effect is not immediately obvious. 

RBC occurs when the temperature difference AT = 
Tb — T t between the bottom and top plates exceeds a 
critical value. At much larger AT the fluid flow becomes 
turbulent. Here we report results obtained with a cylin- 
drical sample of aspect ratio r = D / L ~ 1.00 (D is the 
diameter and L the height) as a function of the Rayleigh 
number 

agATL 3 

Ra = 1 

nv 

over the range 5 x 10 9 < Ra < 10 11 (a is the isobaric 
thermal expansion coefficient, g the acceleration of grav- 
ity, k the thermal diffusivity, and v the kinematic vis- 
cosity). We used de- ionized water as the fluid, with the 
addition of very small amounts of PEO of nominal molec- 
ular weight 4 x 10 6 g/mole, |l4j at a mean temperature 
of 40°C where the Prandtl number Pr = v/ k of the pure 
fluid is 4.38. 

We focused on the effect of the polymer additive on the 
convective heat transport, as expressed by the Nusselt 
number 

Nu = (2 ) 

AATX y ' 

where Q is the heat current passing through the sample, 
A the sample cross-sectional area, and A the thermal con- 
ductivity. As the polymer weight fraction c was increased 
from zero at constant Ra ~ 7 x 10 10 , Nu decreased ap- 
proximately linearly at a rate [1— Nu(c)/Nu(0)]/c ~ 600. 

The measurements were made using the "large sam- 
ple" described in detail elsewhere. [Ill] It had a diameter 
D = 49.7cm and height L = 50.61cm, a cylindrical plex- 
iglas side wall of 6 mm thickness, and aluminum top and 
bottom plates. The plates were magnaplated[l5j to pre- 
vent corrosion and to reduce the deposition of polymer 
from the solution. 
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FIG. 1: The Nusselt number Nu as a function of the polymer 
weight fraction c for Ra — 6.74 x 10 10 . The line is a fit of a 
quadratic equation to the data. 



As shown elsewhere, the finite conductivity of the top 
and bottom plates sli ghtl y reduces the heat transport 
by the turbulent fluid Jill. Il2| For the present purpose 
the uncorrected measurements of Nu are reported. This 
seems appropriate since the change due to the polymer 
addition is of interest and since the plate effect has been 
reported in detail before. 

Measurements were made at constant applied temper- 
ature difference with the mean temperature T m fixed at 
40°C. The maximum applied AT was 10°C, correspond- 
ing to a Rayleigh number Ra — 6.7 x 10 10 . We measured 
the shear viscosity rj of the polymer solutions with an 
Ostwald viscometer and found no change due to the ad- 
dition of the polymer. Thus Ra was calculated using the 
properties of pure water at T m . 

In Fig. [TJ we show the Nusselt number as a function of 
the weight fraction of PEO. The solid line is the result of 
fitting a quadratic equation to the data. It is given by 



Nu(c)/Nu(0) = 1 - 611c(l + 2600c) . 



(3) 



The observation that the dimensionless coefficients in this 
equation are of order 10 3 rather than of order unity re- 
flects the extraordinary sensitivity of this system to the 
polymer additive. 

In Fig. [2] we show all of our measurements in the form 
of the reduced Nusselt number Nu/Ra 03 as a function 
of Ra (using Nu/Ra 03 rather than Nu itself has the 
advantage of being roughly constant and thus of allowing 
greater resolution). The open circles are the results for 
the pure fluid. The solid circles at Ra — 6.74 x 10 10 are 
the same data as those shown in Fig. [TJ The solid squares 
are the results for c = 1.2 x 10~ 4 as a function of Ra. 
One sees that the reduction of Nu due to the polymer- 
decreases slightly with decreasing Ra. 

As mentioned in the Introduction, it is surprising at 
first sight that small amounts of polymer additive have 
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FIG. 2: The reduced Nusselt number Nu/Ra 03 as a function 
of the Rayleigh number on a logarithmic scale. Open circles: 
pure water. Solid circles: Ra = 6.74 x 10 10 and varying poly- 
mer concentration (see Fig. [TJ for a different representation of 
the same data) . Solid squares: 120 ppm polymer as a function 
of Ra. 



any influence at all on RBC. In order to gain some in- 
sight into what may be happening, one needs to consider 
the major features of this system. The temperature dif- 
ference is sustained mostly by two laminar, albeit fluc- 
tuating, thermal boundary layers (BLs), one below the 
top and the other above the bottom plate, with the inte- 
rior of the system nearly isothermal in the time average 
but with vigorous temperature and velocity fluctuations. 
When the BLs are averaged over time and in the horizon- 
tal plane, their thickness is given by A T = X/(2JVu).0 
In our experiment, with L ~ 500 mm, Nu varied ap- 
proximately from 100 to 250 and thus At varied from 
about 2.5 to 1 mm as Ra changed from 5 x 10 9 to 10 11 . 
A useful cartoon of the system is one in which, roughly 
speaking, the BLs adjust their thickness so that they are 
marginally stable, i.e. their Rayleigh number, based on 
Eq. QJbut using a temperature difference equal to AT/2 
instead of AT and a length At instead of L, approxi- 
mately has its critical value for the onset of convection. 
Under these conditions one would expect excitations of 
local flow fields to occur in the BLs at irregular tempo- 
ral and spatial intervals as perturbations act upon the 
BLs. These excitations are emitted from the BLs and 
have become known as "plumes" . It seems likely that the 
polymer additive reduces the perturbations acting upon 
the BLs and thereby reduces plume emission. This would 
result in a reduction of the convected heat transport by 
the system. 

Although both experiment and theory have been lack- 
ing heretofore, very recently^ a theoretical investiga- 
tion, using both direct numerical simulation and ana- 
lytic model calculations, was performed for a RBC sys- 
tem with periodic boundary conditions along the top and 
bottom of the sample. In this system, which seems dif- 
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ficult to achieve experimentally, there are no boundary 
layers and all the dissipation occurs in the bulk. It was 
found that the heat transport is increased significantly by 
the addition of polymers, which is opposite to what we 
found experimentally in the physical system. Since there 
are no rigid boundaries, the theoretical result clearly is 
unrelated to drag reduction in the usual sense. The same 
authors also discuss qualitatively the more realistic case 
of rigid boundary conditions at the top and bottom of 
the sample, and suggest that either an increase or a de- 
crease of the heat transport can be associated with poly- 
mer addition; which case occurs depends on the values of 
relevant parameters. However, a direct comparison with 
experiment does not seem possible at this time. 

In view of the discovery reported in this Letter, a great 
deal more work is warranted on this system. Measure- 
ments of the reduction of Nu as a function of the molec- 
ular weight of the polymer will be illuminating. Optical 
observations of the thermal boundary layers, and of the 
plume density[13[ as a function of polymer concentration, 
will be instructive. Determinations of the Reynolds num- 
ber of the large-scale circulation in this system should be 
revealing. At least for now these investigations unfortu- 
nately are beyond the possible scope of this work. 

This work was supported by the U.S National Science 
Foundation through Grant DMR07-02111. 
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